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Abstract 

In this article, we extend our previous work to study the S-type heavy baryons E c and 
in the nuclear matter using the QCD sum rules, and obtain three coupled QCD sum rules for 
the masses M£ , vector self-energies E„ and pole residues AJ Q in the nuclear matter. Then 
we take into account the effects of the unequal pole residues from different spinor structures, 
and normalize the masses from the QCD sum rules in the vacuum to the experimental data, 
and obtain the mass-shifts 8Ms c = — 123 Me V and 8Ms b ~ — 375 Me V in the nuclear matter. 
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1 Introduction 

The QCD sum rules is a powerful theoretical tool in studying the ground state hadrons both in 
the vacuum and in the nuclear matter, and has given many successful descriptions of the hadron 
properties, such as the masses, decay constants, form- factors, hadron coupling constants, etc [JJ. 
The properties of the light-flavor mesons and baryons in the nuclear matter have been studied 
extensively with the QCD sum rules [21 [21 E] > while the existing works on the heavy quarkonia 
and heavy baryons in the nuclear matter focus on the J/ip, ij c , D, B, D , B Q , D*, D\, A c and A& 
[SI IH1 [7] . The in- medium mass modifications of the Qq and Qqq hadrons differ greatly from the 
corresponding ones of the qq and qqq hadrons due to the appearance of the heavy quark, the full 
propagators of the heavy quarks in the nuclear matter undergo much slight modifications compared 
with that of the light quarks. The upcoming FAIR (facility for antiproton and ion research) project 
at GSI (heavy ion research lab) provides the opportunity to extend the experimental studies of the 
in-medium properties of the mesons and baryons into the charm sector [8l [9J , we can study the 
heavy hadrons in nuclear matter with the QCD sum rules and make predictions to be confronted 
with the experimental data of the CBM (compressed baryonic matter) and PANDA collaborations 
in the future. 

The scattering amplitude of one-gluon exchange can be rephrased into an antisymmetric an- 
titriplet 3 C and an symmetric sextet 6 C in the color-space. The attractive interaction in the 
antisymmetric antitriplet favors the formation of the diquark states in the color antitriplet 3 C , the 
most stable diquark states, the A-type diquark states, maybe exist in the color antitriplet 3 C , flavor 
antitriplet 3/ and spin singlet l s channels due to Fermi-Dirac statistics [TO], and the next stable 
diquark states, the S-type diquark states, maybe exist in the color antitriplet 3 C , flavor sextet 6/ 
and spin triplet 3 S channels jllj . In the heavy quark limit, the heavy baryons can be classified 
as the A-type or S-type baryons according to the spin structures of the two light quarks [T2] . 
In Ref.[7], we study the A-type heavy baryons A c and A& in the nuclear matter using the QCD 
sum rules, and obtain the in-medium positive mass-shifts SMa c — 51MeV and SM\ b = 60MeV, 
respectively. In this article, we extend our previous work to study the properties of the S-type 
heavy baryons E c and in the nuclear matter using the QCD sum rules. 

The article is arranged as follows: we study the heavy baryons S c and in the nuclear matter 
with the QCD sum rules in Sec. 2; in Sec. 3, we present the numerical results and discussions; and 
Sec. 4 is reserved for our conclusions. 
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2 The in- medium Eq baryons with QCD sum rules 



We study the E c and Ef, baryons in the nuclear matter with the following two-point correlation 
functions 11(g), 

11(g) = ifd*xe^ x (HF \T{j{x)J(0)}\iSf ), 

J(x) = e^ujWC^dj (x)^ l5 Q k (x) , (1) 

where the i, j, k are color indexes, Q — c,b, the C is the charge conjunction matrix, and the |\&o) 
is the nuclear matter ground state. The correlation functions 11(g) can be decomposed as 

11(g) = n s (g 2 ,g- U ) + n 9 (g 2 ,g- M )^+n it (g 2 ,g- M )^, (2) 

according to Lorentz covariance, parity and time reversal invariance [2J[3]. In the limit = (1, 0), 
the component IT(g 2 , q ■ u) reduces to IT(go, q), where i = s,q, u. 

We insert a complete set of intermediate heavy baryon states with the same quantum numbers as 
the current operators J(x) into the correlation functions II (p) to obtain the hadronic representation 
PQ, then isolate the ground state Eq baryon contributions, and resort to the dispersion relation to 
rephrase the three components of the correlation functions IT(go, q) in the following form: 

u lM = -i r tjwtb*, (3 ) 

2m J_ 00 u-qo 



where 



A* 2 M* 

AII>, q) = -2ni [8(u - E q ) - S(uj - E q )} , 

A* 2 

AIL>, q) = -27ri-g [6(w - E q ) - 5(uj - E q )\ , 
A* 2 E 

AU u (uj, q) = +2™-|§^ [6(w - E q ) - 8{u - E q )} , (4) 



E* = yM£ 2 Q + g 2 , E q = T, v + E*, E q = E„ - E*, the M£ Q , Su and Af, Q are the masses, vector 

self-energies and pole residues of the Eq baryons respectively in the nuclear matter. 

We carry out the operator product expansion in the nuclear matter at the large space-like region 
g 2 <C 0, and obtain the spectral densities at the level of quark-gluon degrees of freedom, then take 
the limit = (1, 0), and write the three components IT(go, q) in the following form [2J [3J: 

Mq^q) = ^Q(g ,g)(O„) pjv , (5) 

n 

where the C^(go,g) are the Wilson coefficients, the (O n ) PN , which are defined as (\E'o|0„| , I'o)i 
are the condensates in the nuclear matter and can be decomposed as (O) + pn(0)n at the low 
nuclear density, the (O) and (0)n denote the vacuum condensates and the nuclear matter induced 
condensates, respectively. The imaginary parts of the QCD spectral densities can be obtained 
through the formula AIT(w, q) — limit e ^o [IT(w + ie, q) — Tli(u) — ie, q)]. 

Finally, we match the hadronic spectral densities with the QCD spectral densities, and multiply 

both sides with the weight function (oj — E q )e ~~m i , perform the integral J_° duj, 

dwAIIj(w, q)(oj — E q )e~~M I , (6) 
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to exclude the negative-energy pole contributions, and obtain the following three QCD sum rules: 
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where Eq = + if, Eq — m 2 Q + q 2 , xi — 7 ^f-, Sg = uj 2 = s - <f, the w is the threshold 

parameter, and Xi —> in the spectral densities where the function S(s — Eq) appears. We can 
obtain the masses M£ Q , vector self-energies and pole residues A Sq in the nuclear matter by 
solving the three equations with simultaneous iterations. 

3 Numerical results and discussions 

The input parameters of the QCD sum rules in the nuclear matter are taken as (q'q) PN — §/?iv, 
(qq) PN = (qq) + T^f^PN, (««) = (-0.23 GeV) 3 , m u +m d = 12MeV, a N = 45 MeV, (^) PN = 
(^) -0.65GeV PJV , (fi^) = (0.33 GeV) 4 , (qUD iD q) PN + ±(q\g s cxGq) PN = 0.031 GeV 2 p N , 
{qiD iD q) PN + \{qg s oGq) PN = 0.3GeV 2 p N , (qg s aGq) PN = {qg s crGq)+3.0GcV 2 p N , (q^g s aGq) PN = 
-0.33 GeV 2 PN , {qUDoq) PN = 0.18GeVp N , (qg s aGq) = m 2 (qq), m 2 = 0.8 GeV 2 , (qq) 2 pN = 
f(qq) PN x (qq) PN + (1 - f)(qq) x (qq), f = 0.5 ± 0.5, <f = (0.27 GeV) 2 , and p N = (0.11 GeV) 3 

Hi- 

We recover the QCD sum rules in the vacuum by taking the limit p^ = 0, then differentiate 
the Eqs.(7-8) with respect to respectively, and eliminate the pole residues A Eq (here we smear 
the star * to denote the pole residues in the vacuum), then obtain two QCD sum rules for the 
masses M Eq , one can consult Ref.[13j for the technical details. 

In the conventional QCD sum rules pQ, there are two criteria (pole dominance and convergence 
of the operator product expansion) for choosing the Borel parameter M 2 and threshold parameter 
so- We impose the two criteria on the heavy baryon states Eq, and choose the threshold parameters 
s = (10.0 ± 0.5) GeV and (43.5 ± 0.5) GeV, the Borel parameters M 2 = (1.9 - 2.7) GeV and 
(4.8 — 5.6) GeV for the heavy baryons E c and respectively [13]. Finally we obtain the values 
of the masses and pole resides M Ec = (2.54 ± 0.15) GeV and (2.42 ± 0.20) GeV, M E(1 = (5.96 ± 
0.10) GeV and (5.73 ± 0.16) GeV, A Ec = (5.4 ± 1.4) x 10~ 2 GeV 3 and (3.6 ± 1.0) x 10~ 2 GeV 3 , 
A Eb = (8.0 ± 1.7) x 10~ 2 GeV 3 and (5.0 ± 1.2) x 10~ 2 GeV 3 from the QCD sum rules with respect 
to the spinor structures 4 and 1, respectively |13) . 

The experimental values of the masses are M E ++ = (2454.03 ± 0.18) MeV, M E+ = (2452.9 ± 
0.4) MeV, M s o = (2453.76±0.18)MeV, M E+ = (5807.8 ±2.7) MeV and M s - = (5815. 2±2.0) MeV 
from the Particle Data Group [H], the average values are M Ec = 2.454 GeV and M Et = 5.812 GeV, 
respectively, which are consistent with the new experimental data from the CDF collaboration [15] . 
The predicted masses from both spinor structures fa and 1 can reproduce the experimental data 
approximatively, however, the pole residues from the spinor structures /[ and 1 differ from each 
other greatly, i.e. the values have the hierarchy A E( ^ 3> A Eq , where the upper indexes denote 
the spinor structures. We can draw the conclusion that the two QCD sum rules in Eqs.(7-8) in 
the limit pn — can be satisfied by the approximate equal masses but unequal pole residues. If 
we obtain the masses M Eq by dividing Eq.(8) with Eq.(7) in the limit pm = 0, the predictions 
M Ec = 1.40^0 05 GeV and M Ei) = 3.56^q ' X q GeV are much smaller than the experimental data due 
to the unequal pole residues. We can multiply the smaller masses by some coefficients to offset the 
effects of the unequal pole residues and reproduce the experimental data. 

Take the same Borel parameters and threshold parameters as the QCD sum rules in the vac- 
uum [13] , we can obtain the hadronic parameters M E = 1.331q'q3 GeV, -^s = 3.33^q'q7 GeV, 
A Sc = 2.46±°; 22 x 10- 2 GeV 3 , A Efc = 1.25±°;°® x 10- 2 GcV 3 , Ej« = 0.446+°;° 3 7 GeV, E=> = 
0.776t^35 GeV from the three coupled QCD sum rules in Eqs.(7-9). In the limit p^ = 0, we 
can obtain the values M Ec = 1.40±°;^GeV, M Eb = 3.56^* GeV, A Ec = l.99t° f 6 x 10~ 2 GeV 3 , 
A Eb = 8.73lg'g5 x 10~ 3 GeV 3 . In calculations, we have taken the assumption that the pole residues 
A Eq in the QCD sum rules (see Eqs.(7-9)) have the same values. In fact, those QCD sum rules can 

be satisfied with approximately the same in- medium masses M Eq and vector self-energies S„ Q , 
but different pole residues A s , which can be denoted as A E (1), A E (2) and A E (3) from the 
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QCD sum rules in Eqs.(7-9) respectively. We draw this conclusion tentatively based on the QCD 
sum rules for the heavy baryon states in the vacuum [13], and normalize the masses from the QCD 
sum rules in the vacuum to the experimental data to study the mass modifications in the nuclear 
matter. 



"'v 



*Mx Q = — ^ -xM- p , (10) 



and obtain the central values 8My, c = —123 MeV and 8Mj> b — —375 MeV, the SMs Q are the scalar 

self-energies Q . The Eq baryons have a heavy quark besides two light quarks, the heavy quark in- 
teracts with the nuclear matter through the exchange of the intermediate gluons, the contributions 
from the gluon condensates are of minor importance and the modifications of the gluon conden- 
sates in the nuclear matter are mild, we expect the ratios of the mass-shifts SA1^ Q /M^ Q are smaller 
than that of the nucleons. From Fig.l, we can obtain the mass differences AMs Q = Mj Q — Ms Q , 
AM Sc = -0.07 ± 0.02 GeV, AM St = -0.24 ± 0.04 GeV, the uncertainties of the mass differ- 
ences originate from the Borel parameters are about 29% and 17%, respectively. The ratios are 
A , A f Sc = —5.0% and A , A ? Sb — —6.7%, the mass modifications are rather small. 

i"E t M ^b 

If we take into account the uncertainties of the heavy quark masses and threshold parameters, 
Sm c = ±0.1 GeV, 5m b = ±0.1 GeV, 5s^ c = ±0.5 GeV 2 , 5a% b = ±0.5 GeV 2 [13], the values of the 
mass differences AMs Q survive approximately, no additional uncertainties are introduced. At the 
interval / = ~ 1, the masses M|i decrease monotonously with the increase of the parameter /, 
the uncertainty Sf — ±0.5 leads to the uncertainties ±0.06 GeV and ±0.21 GeV for the masses M|j 
and M£ respectively in the nuclear matter, and the corresponding uncertainties for the mass-shifts 
5M Sc and SM Sb are ±105 MeV and ±343 MeV, respectively. 

If we take the Ioffe current to interpolate the proton, the QCD sum rules indicate that there 
exists a positive vector self-energy E^ = (0.23 — 0.35) GeV with the typical values of the rele- 
vant condensates and other input parameters, which is consistent with the values of the vector 
self-energies in the relativistic nuclear physics phenomenology, on the other hand, although the 
scalar self-energy depends strongly on the in-medium four-quark condensate and the nucleon a 
term, a reasonable negative scalar self-energy can be obtained with the suitable parameters [3]. 
There exists significant cancelation between the scalar and vector self-energies, which leads to a 
quasinucleon energy close to the free-space nucleon mass and satisfies the empirical expectation 
that the quasinucleon energy is shifted only slightly in nuclear matter relative to the free-space 
mass. The in-medium self-energies E^ and E^ can be written as 

v at _ 8it 2 <jnPn 
2j — 



M 2 (m u + m d ) ' 
32n 2 p N 



M 2 



(11) 



in the leading order approximation, E^/E^ re —1, which indicates a substantial cancelation 
between self-energies E^ and E^ in the nuclear matter. The mean-field models predicate that the 
typical self-energies of the nucleons in nuclear matter saturation density are E^ re —350 MeV and 
E^ re +300 MeV respectively, which correspond to the real energy-independent optical potentials 
S and V, and significant cancelation between the potentials occurs, the effective non-relativistic 
central potential is about tens of MeV. If the same mechanism works for the E c baryon, the 
vector self-energy should be E^ c re +123 MeV rather than +446 MeV according to the unequal 
pole residues (1) ^ A^ Q (2) ^ (3), the total self-energy E^ c + E^ c re under the condition 
^Sq(3)/As q (2) ~ 1-9' then the quasi-E c energy in the nuclear matter close to the free-space 
E c mass. The present prediction of the mass-shift 8My, c = —123 MeV can be confronted with 
the experimental data from the CBM and PANDA collaborations in the future [8] |9] , where the 
properties of the charmed baryons in the nuclear matter will be studied. 
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Figure 1: The masses from the QCD sum rules in the vacuum and in the nuclear matter versus 
the Borel parameter M 2 , the (a) and (b) denote the E c and E{, baryons, respectively. 



In the non-relativistic harmonic-oscillator potential model, the spectrum of the bound states 
(the energies E n and the wave- functions ip n (x)) and the exact correlation functions are known 
precisely. In Ref.[16], Lucha, Melikhov and Simula try to fit the effective threshold parameter so 
as to reproduce both the ground energy Eq and the pole residue Rq (■0o(O)'0o(O)), or reproduce the 
ground energy E only and take the pole residue i?o as a calculated parameter. They observe that 
the pole residue Rq is determined to a great extent by the continuum contributions, and draw the 
conclusion that the ground-state parameters extracted from the QCD sum rules have uncontrolled 
systematic errors if the continuum contributions are not known exactly and modeled by an effective 
continuum. In the real QCD world, the hadronic spectral densities are not known exactly, the pole 
residues or decay constants in some (or most) cases are not experimentally measurable quantities, 
and should be calculated by some theoretical approaches, the true values are difficult to obtain. 
On the other hand, the hadronic spectrum densities in the non-relativistic harmonic-oscillator 
potential model are of the Dirac 6 function type even in the limit n — > oo, while in the case of the 
QCD, the widths of the higher radial excited states become broader gradually before submerging 
into the asymptotic quarks and gluons. For example, the widths of the tt, 7t(1300), 7r(1800), • • • are 
- GeV, (0.2 - 0.6) GeV, 0.208 ± 0.012 GeV, • • • respectively, and the widths of the K 7 #(1460), 
#(1830), ••• are - OGeV, - (0.25 - 0.26) GeV, - 0.25 GeV, ••• respectively [Hj- We cannot 
estimate the unknown systematic uncertainties of the QCD sum rules before the spectral densities 
in both the QCD and phcnomcnological sides are known with great accuracy. In the present case, 
the situation is even worse, the theoretical works on the Eq baryons in the nuclear matter are rare, 
and no experimental data exist. 



4 Conclusion 

In this article, we extend our previous work on the A-type heavy baryons Aq to study the prop- 
erties of the E-type heavy baryons Eq in the nuclear matter using the QCD sum rules, and 
derive three coupled QCD sum rules for the masses M£ Q , vector self-energies E„ and pole residues 
A|; in the nuclear matter, then obtain the values M£ c = 1.331q q3 GeV, M£ b = 3. 33^07 GeV, 
A^ c = 2.46±g;?l x 10- 2 GeV 3 , A£ b = 1.25±8;j£ x 10- 2 GeV 3 , E^ = 0.446±g;g^ GeV, E^ = 
0.776t^3 2 GeV. In the limit p N = 0, the predictions M Sc = 1.40±g;g| GeV and M s „ = 3.56±£;" GeV 
are much smaller than the experimental data due to the unequal pole residues from different spinor 
structures 7 • q and 1. We normalize the masses Ms c = 1.40 GeV and M^ b = 3.56 GeV from the 
QCD sum rules in the vacuum to the experimental data, and obtain the mass-shifts in the nuclear 
matter <5Me c = —123 MeV and SMs b = —375 MeV, which can be confronted with the experimental 
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data in the future. 
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